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Abstract 

We study the neutrino oscillation problem in the frame- 
work of the wave packet formalism. The neutrino state is de- 
scribed by a packet located initially in a region S (source) and 
detected in another region D at a distance R from 5*. We ex- 
amine how the oscillation probability as a function of variable 
R can be derived from he oscillation probability as a func- 
tion of time t, the latter being found by using the Schrodinger 
equation. We justify the known prescription t ^ R/c without 
referring to a specific form of the neutrino wave packet and 
only assuming the finiteness of its support. The effect of the 
oscillation damping at large R is revealed. For an illustration, 
an explicit expression for the damping factor is obtained using 
Gaussian packet. 
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1 Introduction 



While in the minimal formulation of the standard electroweak model 
neutrinos are massless, there have been strong experimental evidences 
for nonzero neutrino masses and mixing of neutrino flavors which for- 
tify the idea of Pontecorvo [1] that nonzero neutrino mass and lepton 
flavor violation would produce spontaneous oscillation of neutrinos 
from one flavor to another in a manner similar to the strangeness 
oscillations of neutral kaons. The consistent quantum-mechanical 
treatment of neutrino oscillations is important for interpretation of 
the current and future experiments. 

A typical neutrino oscillation experiment can be described by the 
following scheme. There is a neutrino source which produce neutrinos 
of a definite fiavor Va (a = e, /x, r) in a finite region S. The detector 
of neutrinos of the same or another flavor vp is located in another 
finite region D at a distance R from S. The regions S and D are 
supposed to be well separated from one another. The probability to 
detect with /3 7^ a is in general nonzero. The theoretical task is 
to predict the R dependence of this probability. For simplicity, only 
the neutrino oscillations in vacuum are considered here. 

There are many approaches to this problem. The simplest is to 
find first the temporal evolution of the neutrino state and then to 
convert it into the spatial evolution. Usually it is assumed that the 
neutrino state at t = has a definite momentum that is described 
by a plane wave. Then, by applying the Schrodinger equation, one 
finds the probability W{t) to detect the same neutrino fiavor at i > 
(survival probability) or another flavor (transition probability) aXt > 
(Sect. [5]). In order to convert the time dependence of W into the R 
dependence, one apply the classical propagation prescription (in the 
terminology of Ref. [2]), according to which t is simply replaced by 
R/c (where c = 1 is the velocity of light), with R being displacement 
of the classical ultrarelativistic particle at time t. The deficiencies 
of this prescription are well known (see, e.g., Ref. [3, Ch. 2.3.3]). In 
particular, the formula R = vt is only valid for a classical particle 
trajectory and is inconsistent with quantum postulates. 

It is generally accepted that one must describe the neutrino state 
by a wave packet rather than by a plane wave. The relevant consid- 
erations are called in Ref. [2] intermediate wave packet models. Many 
aspects of of these models were discussed in Refs. [4-17]. The present 
paper belongs to this class. Some peculiarities of our approach are 
presented in the concluding section (Sect. [4]). One of our aims (in a 
sense, pedagogical) is to propose an example of how one can obtain 
the time-to-space conversion within the framework of the quantum 
postulates. 

In other formalisms, like the external wave packet models and 
quantum field theoretical or S matrix approach (see Ref. [2] and ref- 
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erences therein) , the wave packet description is used for the particles 
which either create neutrino (e.g., a pion, when neutrinos are pro- 
duced through the 7r^2 decay) or detect it (e.g., an oxygen nucleus in 
a water-Cherenkov detector). Being more advanced, these approaches 
are more cumbersome and usually are aimed at experiments of a spe- 
cific type. The intermediate wave packet models are simpler and, in 
a sense, more general, being applicable for generic neutrino source 
and detector. What is more, they are sufficient for the purpose of 
the present study: to show how the flavor oscillations in time leads 
to these in space (Sect. [31). 

Our main results are summarized in Sect. |3) 

2 Neutrino oscillations in time 

In order to define what "neutrino oscillations in time" means and to 
prepare equations necessary for the next Section, let us discuss the fol- 
lowing problem (c.f., for example, Ref. [3]). One free electron neutrino 
state with definite momentum k {weak eigenstate) is prepared at the 
moment t = 0. Let \l/(0) = |e,k) denotes this initial state (as usual, 
we omit spin indices). The evolution of the state in time is governed 
by a time-independent and nondiagonal in the flavor basis Hamilto- 
nian H, the propagation Hamiltonian, according to Ref. [2, Ch. 2]. 

Then ^'(i) = cxp (^—iHtj ^'(0). According to the Pontecorvo's sup- 
position [1], H has eigenvectors |j, k) {j = 1,2,...) with definite 
value of the total momentum k (mass eigenstates). The correspond- 
ing eigenvalues are equal to rUj at k = and to Ej = y'jkp -I- m| at 

k 7^ 0. The weak eigenstate VE'(O) may be expanded in terms of the 
mass eigenstates (below, we assume j = 1, 2 for simplicity): 

|e,k) = cos6'|I,k) +sin6'|2,k). (1) 

Then 

^{t) = \e,k;t) ^ e-*^i*cos6'|I,k) +e~'^^*sin6l|2,k). (2) 

This equation provides all information on the neutrino state at time 
t (expectation values of observables, probabilities). In particular, 
the probability of finding the system in the initial state (survival 
probability) is given by 

W^k(i) = I (e, k|e, k; t) I' = 1 - i sin^ 29 [1 - cos(S2 - Ei)t] . (3) 

For ultrarelativistic neutrinos (|kp ^ 2)? Eq. ([3]) is simplified to 



2'Kt 



(4) 
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with 

y _ 47r|k| 

\ml -ml\' 

This probabihty is a periodic function of time (with the period T) 
and provides the simplest example of neutrino oscillations in time. 
Replacing thy R/cva Eq. ([4]) one obtains the "standard" formula for 
the ultrarelativistic neutrino oscillations in space: 

W^k(i?/c) ~ 1 - isin^ 261 

with 

A = Tc. 

Note 1. For deriving Eq. ^ we imply the unit normalization of 
the mass eigenstates: 

(i,k|j,k) = 5^J. 

Since these states are the eigenvectors of the momentum operator P, 
this assumption means that the operator P has a discrete spectrum. 
This assumes that the system under consideration is in a large spa- 
tial volume and the usual periodicity conditions are imposed (or the 
opposite boundaries of the volume are identified). However, in what 
follows we will approximate the summations over k by integrations 
over the 3-momentum. 

Note 2. Let us consider the muon neutrino created in the decay 
TT"*" /i+ + t'p, 7r+ being at rest. From the energy- momentum con- 
servation it follows that the neutrino state is a superposition of the 
plane waves jkj), j = 1,2, ... , with the |m unambiguously defined 
by the masses m^r, and rrij (see, e.g., Ref. [19, Sect. 3]). In con- 
trast, the "Pontecorvo's prescription" ^ assumes that the neutrino 
initial state is a single plane wave with an arbitrary momentum. Here 
we accept for simplicity just this prescription. 

3 Neutrino oscillations in space 

In any neutrino oscillation experiment the source of neutrinos with 
definite flavors is localized in a finite volume S and the detector of 
neutrinos (of the same or different flavor) is localized in another finite 
region D placed at a distance R from S. The source cannot prepare 
a neutrino with definite momentum k (plane wave). The emitted 
neutrino (e.g., electron neutrino) must be described by a packet, i.e. 



cos 



f2-KR 



(5) 
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by a superposition of the states |e,k), Vk which form a complete set 
of states: 

|e,5) = J d3fcg(k)|e,k), I d'k\~g{k)\' = 1. (6) 

In the coordinate representation, the state \e,g) is described by the 
fimction 

5(x) = (e, x|e, g) = (2^)-3/2 J d^kg{k) exp (ikx) , (7) 

which is the amphtude of the density probabihty that electron neu- 
trino in the state \e,g) has the position x. It is assumed that (7(x) 
is localized in S in the sense that g(x) is negligibly small (physically 
indistinguishable from zero) outside S. 

After multiplying both sides of Eq. ([T]) by g(k) and integrating 
over k we obtain the expansion of \e,g) over the mass eigenstates: 

\e,g) ^cose J d^kg{k)\l,k) +sme J d^kg{k)\2,k) . 

Analogously, from Eq. ^ one gets 
\e,g;t) = exp {-iHtj \e,g) 

= cos0 J kg {k)e-'^'' \1, k) + sin J d^kg{k)e-''^^'\2,k), 

where 



1 2 I rtn 2 



= ^|k|- +mj. 

For the amplitude to find neutrino in the point x we obtain 
(e, x|e, g; t) = cos^ %i(x, t) + sin^ 9g2i^, t), 

c/j(x, t) = (27r)-3/2 J (fkg{k) exp [^i {Ejt - kx)] . (8) 

Here we have used the equalities 

(e,x|l,k) =cos6le'''^, (e,x|2,k) = sin6le^''''. 

The functions gi (x, t) and g2 (x, t) describe the evolving wave pack- 
ets of particles with masses mi and m2, respectively. According to 
Eqs. ([6]) and ([8]), they coincide with the packet ^(x), Eq. ([7]), at the 
initial moment t = 0. It is shown in Appendix A that, neglecting the 
spread of the packet, the evolution of .gj(x, t) reduces, up to a phase 
factor (close to 1 in the ultrarelativistic case), to the "drift" of the 
initial packet g(x) with the velocity Vj — k/Ej: 
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In other words, the wave packets gj{'x.,t) approximately follow the 
trajectory of a classical particle. 

Let us now assume that the detector D registers the electron 
neutrinos. The probability to find in the state \e,g;t) one electron 
neutrino localized within the region D is equal to 



(fx\{e,K\e,g;t)y 



= J d^x^^cos'^0\gl{x,t)\^ +sm'^9\g2{x,t)f' 
+2 cos^ e sin^ e Re [gl (x, t)g2 (x, t)]} . 



(9) 



This probability may also be represented through the average value 
of the projection operator 



Td = d a::|e,x)(e,x|, 



D 



namely 



Wg{t) = 



{e,g;t\fD\e,g;t) 



We consider now the behavior of Wg as a function of t. Let R be 
the vector connecting the center of S with the center of D {R = |R-| 
being the distance between S and D), see Fig. [1] Suppose that the 




Figure 1: The mutual disposition of the neutrino source S and neu- 
trino detector D. The dimensions of the source and detector and the 
distance R between S and D are implied to satisfy the inequalities 
Ls <^ Ld R- The space between dashed lines is covered by the 
supports of moving wave packet as it shifts and (slowly) spreads. The 
vector p is the average momentum of the packet. 

module p of the average momentum of the packet is much greater 
than TOi and TO2 (the neutrino are ultrarelativistic) and therefore 
the group velocities vi and V2 are both approximately equal to the 
velocity of light. 

So, the packets 51 (x, t) and 32 (x, t) enter the region D (and leave 
it) almost simultaneously. Their nonoverlapping will be discussed 
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later on. Till the moment t — R/c, the moving supports of the 
packets gi and 32 do not cover the region D (supp(gj) n D = 0) and 
therefore Wg{t) = if i < R/c. After the moment {R + Ld)/c (where 
Ld is the dimension of the region D in the direction of the vector R) 
the packets get out of D and we have Wg{t) = again (see Fig. 
Let us assume further that 

Ls^Ld^R. (10) 

where L5 is the dimension of 
the region S (which may be 
considered as coinciding with 
the dimension of the initial 
packet) along the vector R 
(see Fig. [1]) . Then we can ne- 
glect all marginal effects, e.g. 
the gradual growing of Wg (t) 
as neutrino packets are en- 
tering the detector. So we 
obtain that Wg{t) is nonzero 
only within the time interval r ~ (R/c, R/c + Lo/c). Inside this 
interval, the probability Wg{t) is approximately equal to Wg{R/c), 
Ld being much smaller than R. Note that if the packet spreading is 
taken into account, one cannot consider arbitrary large values of R 
since the support of the spread packet must not exceed D. 

It is sometimes convenient to deal with the number of particles 
rather than the probabilities. Consider an ensemble of Nq neutrinos 
in the initial state |e,g). Then the total number A'^ of neutrinos Ve 
falling into the detector at all times from t = to 00 is proportional 
to NQWg{R/c). Being independent of time, this number depends 
upon R. Since the neutrinos have a nonzero average momentum p 
then, during the time r there is a nonzero current of neutrinos in D. 
For example, we can consider a detector which uses the Pontecorvo's 
reaction 

i^e + CI Ar + (11) 

for catching the electron neutrinos (see, e.g., Ref. [18, Ch. 2.3.3]). 
The number of the resulting Ar nuclei depends on the number N of 
electron neutrinos falling into the detector, density of the chlorine 
nuclei in D, total cross section of the reaction pT|) . etc. However, all 
these factors but A^ do not depend upon R. We may assume that 
all detectors placed at different distances R from S are identical. It 
is pertinent to remark that, assuming the inequalities (jlOp . we imply 
that each neutrino detected in the chlorine detector had been actually 
emitted by a small piece of the Sun (not by the Sun as a whole) . The 
total number of neutrinos detected in the chlorine experiment must 











% 


>- 



R/c R/c+L^c 

Figure 2: The probability of electron 
neutrino localization in the detector 
region as a function of time. 
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be confronted to the probability Wg{R/c), Eq. summed over all 
these small pieces which make up the Sun. 

Let us remind that the chlorine experiment does not register the 
time when Ar nuclei appear while measures the number of the argon 
nuclei accumulated in the detector during a long period. Thus, in- 
stead of the dependence on time there appears the dependence on R 
(through A'^). If R is fixed then the actual time of appearing of Ar 
nuclei in D is also fixed: t = R/c. According to Beuthe [2, Ch. 2.3.1], 
this is called conversion of the evolution in time into the evolution in 
space. This conversion has been obtained above by using the wave 
packets but their specific form was not employed (only finiteness of 
their supports was used). 

Our next task is to calculate numerically the probability Wg, 
Eq. Ill) in order to ascertain if one can obtain from Wg the standard 
formula ^ for Wk{t). This calculation reduces to the calculation of 
the integrals 

/ d^x\g,{^,t)\\ 3 = 1,2 (12) 
Jd 

and 

/ d\gl{-K,t)g2{^,t) = h2{t). (13) 

JD 

Let us assume that the packets 51 (x, t) and 52 (x, t) are entirely inside 
D at t £ T, T = [R/c, R/c + Ld/c) (see note below). Then the 
integrals (|12p can be replaced by the integrals over the whole three- 
dimensional space. The latter do not depend upon t, both being 
equal to the normalization integral 

jSx\g{^f = l. (14) 

Note. Besides shifting in space, the packets 5i(x, t) also spread (see 
Appendix A). There is a boundary value i?spr of the distance R such 
that the spread packets completely cover the region D when R > 
Rspr- Then the integrals ([T2| are not equal to 1. If i? < _Rspr then 
the packet's supports are in D. One may define i?spi. as follows. Let 
gi{:x.,t) are Gaussian, see Eq. ((23|) or Eq. p^)) . Then its transversal 
spatial dispersions equal to cr J 2 (t) , see Eq. (p4)) , its spacial dimension 
being (Ti,2(i)- The dimension is equal to Ljj (the region D dimension) 
when t ^ tspi ~ Rspr/c = 2Ldijp (Eq. (|24p is used and inequality 
Ld ^ Ls ^ a in presumed). Here a is the spatial dimension of the 
initial packet (;(x, 0) and p is the magnitude of the average momentum 
of the packet. So integrals (fT2| equal to 1 if i? < _Rspi- = iLnapc. 
If, for example, 5 is a macroscopical packet (say cr ~ 1 cm) with 
the average energy _E ~ p of the order of several MeV then _Rspi- ^ 
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The integral /12 defined by Eq. (fT^ is more complicated. It ap- 
proximately equals to 1 when t is small (along with R; remind that 
t £ T with T = {R/c,R/c + Ld/c)) because gi(x, ~ g2{x,t) at 
small t. The integral /12 vanishes when t R/c is sufficiently large 
because gi and 32 cease to overlap (being both in D if i? < i?spr) 
due to different group velocities, vi ^ V2- For such R the probability 
Wg{R/c) does not longer depend on R and does not oscillate. This 
phenomenon does not exist in the standard formula (O- 

For a numerical illustration we calculated in Appendix B the in- 
tegral I12 by using the Gaussian packet 



g(k) (X exp — (p — k)^ cr^ 



(15) 



Here p is the average momentum and is the spatial dispersion of 
the packet. It is assumed that p = |p| ^ mi and that a is macro- 
scopically large, so that we have 



per > 1. 



(16) 



We suppose that R <C -Rspr ~ Locrpc so that the integrals are 
equal to 1. The resulting value of Wg{t) at the moment t — R/c can 
be represented as 



Wg{R/c) = 1~ ^sin^ 20 



I — V cos 



27ri? 



(17) 



where 



and 



V 



exp 



IPI^ 



(18) 



A = Tc = 



47r|p| 



Ato^ = \mf — 'm\\. 



Am^ 

So, the spatial oscillations of Wg are damped by the factor T). The 
damping can be neglected when R is such that I? ~ 1 and it is signif- 
icant when R is such that V < 1. Let us designate the corresponding 
boundary value of R as i?damp and define it as that value of R at 



which V ^ 1/e. Then, as it is seen from Eq. (|T8|) . i?damp ~ Aapc 
and, according to the inequality (fT6|) . -Rdamp ^ A. So the damping 
can be neglected, e.g., for R of the order of one or several oscillation 
length A. In this case, Wg{R/c) turns into the r.h.s. of Eq. (O for 
Wk{R/c), i.e., into the standard neutrino oscillation formula. 

We have 



Now, let us compare -Rdamp with i?spr- 



-Rd 



amp 



Rs 



A 



P 



Am^ 



1 



P 



1 



(19) 
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Due to the inequality p ^ rrii the first factor in r.h.s. of Eq. (|19p 
is large, but LdV Ato^ can be estimated as being much larger and 
therefore i?damp < -Rspr- 

Thus we have obtained the following. When R ^ i?damp the 
standard formula is correct. As R grows and becomes > i?damp, 
the damping turns out to be essential and the standard formula fails, 
while our result, given by Eq. ^ still holds. For very large R > -Rspr, 
our formula becomes also incorrect due to the packet spreading. We 
conclude therefore that there exists an interval of R values for which 
our formula ([9|) (containing the damping) is valid, while the standard 
formula is incorrect. 



4 Conclusion 

The above consideration of neutrino oscillations is of pedagogical 
interest since it does not use the reasons inconsistent with quan- 
tum postulates. In particular, we have presented a purely quantum- 
mechanical justification of the "classical propagation prescription" 
t I — > R/c. Usually, this prescription is grounded on the classical 
reasons inconsistent with quantum mechanics. 

Our justification is based on the observation that fixing the de- 
tector at a distance R from the neutrino source S automatically fixes 
the neutrino sojourn time in the detector volume. The standard re- 
sult given by Eq. (O for the spatial neutrino oscillations follows from 
Eq. ([9|) when the spreading of the neutrino packet is neglected. It 
is shown that the spreading cannot be neglected if R is too large. 
In this case the dimensions of the spread neutrino packets become 
larger than the detector dimension Ld- the packets cover the detec- 
tor. An explicit estimation R > Rgpr is obtained by using a specific 
(Gaussian) form of the packet (the meaning of i?spr is explained in 
the Note after Eq. ([H), Sect.[3l). 

We point out another qualitative reason for the breakdown of the 
standard formula. Consider, for example, the case when dimensions 
of the packets gi(x, t) and (72 (x, t) of neutrinos with masses nii and 
7712 are both less then Ld (both being in D) but their supports do 
not overlap due to different group velocities vi ^ V2 (see text before 
Eq. (fT5|) ). Then the oscillations cease. The explicit example with the 
Gaussian packets shows that in the case of incomplete overlapping 
there arises a damping of oscillations which is significant for large R. 
{R ~ i?damp -Rspr; for i?damp sce the tcxt bcforc Eq. ([TO]) ). The 
resulting formula for the spatial oscillations reduces to the standard 
one for R ^ i?damp when the damping can be neglected. 
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Appendix A: Packet evolution 



Let us estimate the integral (H)). Let the function g(k) in Eq. ([8]) be 
neghgibly smah outside the sphere |k— < Ip with the center p and 
radius Ip in the momentum space. We may choose Ip to be several 
times larger than the dispersion of the momentum distribution ^(k). 
Assume that g(k) = S{p — k), where S'(k) is a spherically symmetric 
function which is small when |kp > 1^. The average momentum 
/ d'^kg{k)k is equal to p. 

Having in mind (ultra) relativistic neutrinos one can safely assume 
that p ^ rrii. Nevertheless, for generality, we leave p arbitrary. The 
integral ([5]) will be estimated for the case Ip E ^ ^/p^ + m^. The 
corresponding spatial dimension of the initial packet (7(x, 0) satis- 
fies the inequality 1/ E. In particular, may be macroscopically 
large; e.g., for neutrino energies of the order of several MeV it can be 
as large as 1 cm. 

Substituting ^(k) ~ S{p — k) into Eq .([5]) and changing the vari- 
able k ^ p — k we arrive at the following result (index j is omitted 
for simplicity) 

5(x, t) = (27r)"3/^ j exp(ikx)S'(p - k) exp [-ityj\k\'^ + m?^ 
exp(ipx 



(27r)3/2 



d3fcS'(k)exp (kx + tV|p-kp 



(20) 



Since 5'(k) does not vanish only when |k| < Ip, we may expand 
a/Ip — kp -I- -m? in powers of k/i?: 



irr ~ E 



pk |kp - (vk)- 

^2 + 2E^ 



P 

^=E- 



(21) 



Retaining only the two leading terms in the square bracket of Eq. ([21 
we obtain 



5(x, t) ~ (27r)-3/2 exp [i(px - Et)] I d^kS{k) exp [-ik (x - tv 



g(x- v<). 



(22) 



Here g(x) = ^(x, 0) is the initial packet. Eq. (|22p means that the 
packet evolution in this approximation reduces to the shift of the 
initial packet; the phase factor e^'*™ is inessential. 

The additional spreading of the packet originates from other terms 
in the square bracket of Eq. (|2ip . To give an idea of their effect on 
the packet evolution we shall take into account the third term in the 
bracket assuming that the momentum p is directed along R that is 
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p = {0,0, vE). Moreover, let us consider the particular case when 
the distribution ^(k) is Gaussian: 



^(k) = ^exp (-Ikpd^) = ^exp [-a'^{kf + + kj)] 



(23) 



Here ^ is a constant and is the spatial dispersion of the packet. 
Then g{x,t) reduces to the product of three integrals: 



Here 



g(x, t) = Aexp [i (px - Et)] hhh 



Vi = V2 = 0, V3 = V 

T it 



Oi = 02 = (J 



, it(l - v^) 



2E' 2E 
Note that aj are complex numbers with positive real parts. According 
to Ref. [20, Ch. 2.5.36.1] the integrals Ij are 



I, =2 



due 



cos K [xj — Vjt) 



Hence we obtain 



|g(x,t)|2 = lAhhhl^ = Const • exp 



■ exp 



{xj - Vjty 

4a, 



where 



al{t)^alit) 



al{t)=a' 



t^l~v^) 



(24) 



and the constant can be determined from the initial condition and 
the normalization condition (|14p for the function g(x). 

The packet spatial dispersions along the axes xi, X2 and x^ are all 
equal to at the initial moment, t = 0, and increase with increasing 
time according to Eq. ([24|) . In the ultrarelativistic case (1 — ^ 1) 
and at large enough t, a\ 2(i) may be much larger than cr|(i). Taking 
into account that, according to our assumption, a is macroscopically 
large (cr 3> 1/p), one can also conclude that increase in the packet 
dimension in all directions is much less than the spatial shift vt of 
the packet as a whole. Some numerical examples when the packet 
spreading can be neglected are presented in Ref. [21, Ch. 3.1]. The 
physical conditions at which the spreading can be neglected in our 
task are examined in Sect. [3l 
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Appendix B: Evaluation of the integral I12 



Let us evaluate the integral ([13]) defined in Sect. [Sj for a practically 
important particular case. Assume i € r and R < Rspr (with i?spr 
defined in Sect. [3]) so that the supports of both gi(x, t) and 52 (x, t) are 
in D. Hence the integration over the region D in can be extended 
to the whole space. Now let us use the same assumption about the 
Fourier transforms gi(k) as in Appendix A. Then (cf. Eq. (PO)) ). 



exp(ipx) 



(27r)3/2 

and we can rewrite Eq. (|13p as 

Il2{t) = 



J (fkS{k) exp -i (kx + 



|p-kP 



X exp < it 



(p - k) - - V (p - k) 



(25) 



The square roots in Eq. ([25]) prevents exact calculation of the integral. 
However it can be estimated for the ultrarelativistic case and taking 
into account that the function ^(k) is only appreciable for |k| ^ |p|. 
These simplifications allow us to use the approximation 



(p — k)^ — m\ 



(p-k)^ 

Substituting Eq. (gl]) into Eq, 
it{m\ ~ 7712) 



to5 



2|P| 
we arrive at 



1 



(kp) 



(26) 



h 



exp 



2|P| 



d-^k 



5(k) 



exp 



it{m\ — rrij) (kp) 



2|P|' 



In particular, for the Gaussian distribution given by Eq. (j23[) and 
considering that 



/l2(0) = j d3x|g(x)|' = l, 



we obtain 

h2{t) = exp 



it{m\-m2) ( \-p\t{m\ ~ rn^) 



2|P| 



exp 



where 



2iTTt 

T 



T — 



47r|p| 



|P|(^ 
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